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Abstract 

Interfacing atomistic-based with continuum-based simulation codes 
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is now required in many multiscale physical and biological systems. 
_ * We present the first results from coupled atomistic-continuum simula- 

tions on 190,000 processors. Platelet aggregation in the patient-specific 
I model of an aneurysm has been modeled using a high-order spectral/hp 

i— ( element Navier-Stokes solver with a stochastic (coarse-grained) Molec- 

ular Dynamics solver based on Dissipative Particle Dynamics (DPD). 
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1 Introduction 
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Low resolution movie: I Video ll 
High resolution movie: Video 2 

The fluid dynamics video (in low resolution Video 1 and in high resolu- 
tion Video 2) shows: 
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1. Dissipative Particle Dynamics (DPD) simulation of healthy 
(red) and malaria-infected (blue) red blood cells (RBC). The 

membrane of each RBC is represented by 500 DPD particles [1]. The 
blood plasma is also represented by DPD particles, shown as spheres 
colored by the local velocity magnitude. The ensemble average velocity 
is computed using window proper orthogonal decomposition, and is 
presented at several slices along the computational domain. 

2. Patient-specific domain of major brain arteries including cir- 
cle of Willis, internal carotid, vertebral and communicating 
arteries; as well as large aneurysm. At the four inlets patient- 
specific inflow boundary conditions have been imposed. At the multi- 
ple outlets the pressure is regulated by the RC boundary conditions [2] . 
The main flow characteristics are Re = 394 and Ws = 3.7. Computa- 
tional domain is subdivided into four overlapping patches composed of 
about 0.5M spectral elements. In each element the solution is approx- 
imated with six-order polynomial expansion. Instantaneous velocity 
field is represented by streamlines, which also highlight the presense 
of an unsteady vortex inside the aneurysm. 

3. Platelet aggregation at the wall of the aneurysm. Simulation 
is performed using a coupled continuum-atomistic solver. The con- 
tinuum solution in the macroscopic domain is obtained by Nektar, 
while the atomistic solution inside 4 mm 3 microscopic subdomain is 
computed by DPD-LAMMPS. The velocity boundary conditions for 
the microscopic domain are updated every 200 time-steps using data 
simultaneousely computed by Nektar. Platelet aggregation simula- 
tions have been performed on up to 190,740 processors of CRAY XT5 
at ORNL and on up to 131,072 processors of BlueGene/P at ANL. 
Platelet aggregation modeled is described in [3]. More information on 
the coupled continuum-atomistic solver will be provided in the SC'll 
Gordon Bell Prize finalist paper [4]. 
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